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Qh! Abstract 
(D . 

' Topological strings on Calabi-Yau manifolds are known to undergo phase transitions 

at small distances. We study this issue in the case of perturbative topological strings 
^ ■ on local Calabi-Yau threefolds given by a bundle over a two-sphere. This theory can be 
regarded as a q-deformation of Hurwitz theory, and it has a conjectural nonperturbative 
description in terms of q-deformed 2d Yang-Mills theory. We solve the planar model and 
find a phase transition at small radius in the universality class of 2d gravity. We give 
strong evidence that there is a double-scaled theory at the critical point whose all genus 
free energy is governed by the Painleve I equation. We compare the critical behavior 
of the perturbative theory to the critical behavior of its nonperturbative description, 
which belongs to the universality class of 2d supergravity, and we comment on possible 
implications for nonperturbative 2d gravity. We also give evidence for a new open/closed 
duality relating these Calabi-Yau backgrounds to open strings with framing. 



^Also at Departamento de Matematica, 1ST, Lisboa, Portugal. 
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1 Introduction 



It has been recognized for a long time that string theory provides in a natural way a 
deformation of classical geometry. The deformation parameter is given by 



where tg is the lenght of the string and R is the characteristic size of the target space. 
When this parameter is very small, the target geometry can be regarded as a classical 
background corrected by stringy effects. However, when this parameter is big, classical 
geometric intuition breaks down and one is forced to use some notion of stringy or quantum 
geometry. 

These phenomena have been studied in detail in the context of type A topological 
string theory on Calabi-Yau manifolds, where the characteristic size R of the target 
space is set by the Kahler moduli (see ,2] for reviews and a list of references). When 
these sizes are small, so that the deformation parameter (11. ip is big, the breakdown 
of classical geometry can be made precise quantitatively by looking at the behavior of 
topological string amplitudes. Typically, the large radius expansion of the amplitudes 
becomes divergent at a critical value of the Kahler moduli, and this divergence signals the 
onset of the stringy regime. We will refer to this change of regime as a phase transition. 
In known cases (like the quintic), this transition exhibits a universal behavior which has 
been identified with the c = 1 string at the self-dual radius [3l H] . 

In this paper we will analyze the critical behavior of type A topological string theory 
in the local Calabi-Yau manifold 

= 0{v - 2) © 0(-p) (1.2) 

This theory has been studied from a mathematical point of view in [5] and it is closely 
related to Hurwitz theory. From a physical point of view, it was proposed in [6j that the 
model could be defined non-perturbatively by using a q-deformed version of 2d Yang- 
Mills theory. Although perturbative topological string theory on the background 
shares many known properties with other local Calabi-Yau manifolds (like for example 
Gopakumar-Vafa integrality), we will show in this paper that for p > 2 it exhibits a phase 
transition at small radius which is not in the same universality class as the examples 
considered so far, i.e. the c = 1 string. Rather, the critical behavior is in the universality 
class of pure two-dimensional gravity (i.e. the (2, 3) model) for all p > 2. 

It turns out that, as in the case of matrix models, the critical behavior of higher genus 
amplitudes makes it possible to define a double-scaled theory at the transition point. 
The double-scaling limit erases a lot of information about the original topological string 
theory, since it only keeps the leading singularity in the genus g free energies, but it is 
interesting for at least three reasons: 
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(1) It gives a compact characterization of the critical behavior at every order in string 
perturbation theory. 

(2) It makes possible to study the free energy at all genus, and provides information 
about the nonperturbative structure of the theory. 

(3) As argued in [TIE] in the context of unitary matrix models describing thermal AdS, 
the phase transition might be an artifact of perturbation theory, and the double-scaling 
limit provides a way to smooth out the critical behavior. 

We give convincing evidence that the double-scaled free energy is described for all p > 
2 by the Painleve I equation which characterizes 2d gravity [Hllinilll]- The resulting scaled 
free energy is known to be non-Borel summable as a genus expansion, and this confirms 
that perturbative topological string theory on these backgrounds lacks fundamental non- 
perturbative input. Interestingly, the non-perturbative definition of the model proposed 
in [6] also exhibits critical behavior [I2l |T3[ [H] , but this time in the universality class of 
the Gross- Witten-Wadia unitary model [T6] described by Painleve II [171 [18] , which 
is also the universality class of 2d supergravity pTO]. The Painleve II equation is known 
to have a well-defined unique solution with the appropriate boundary conditions |2Uj . 
and this can be regarded as further support for the proposal of [6] as a non-perturbative 
definition of the theory. 

Conversely, our result shows that 2d gravity and 2d supergravity can be respectively 
embdedded in a one-family of topological string theories and their nonperturbative holo- 
graphic completions. The relation that emerges between 2d gravity and 2d supergravity 
can be better understood by using the results of [611 165], and we speculate that the pro- 
posal of [6] , together with our results on the critical behavior, give further support to the 
nonperturbative solution of 2d gravity presented in |64j . 

Mathematically, perturbative string theory on Xp can be regarded as a q-deformation 
of Hurwitz theory, and the Gromov-Witten invariants of Xp promote simple Hurwitz 
numbers to polynomials. This ingredient is very important in our study, and indeed the 
critical behavior we describe for p > 2 is already present in the Hurwitz case. This was 
noted at the planar level in [21], and it is implicit in the results of [22] for higher genus 
Hurwitz theory. Our results are also probably related to the description of topological 
gravity in terms of the asymptotics of Hurwitz numbers due to Okounkov and Pandhari- 
pande [23]. In fact, we can reach higher critical points by perturbing the model with 
higher Casimir operators. This perturbed model is an analog of the general one-matrix 
model, while the theory of Hurwitz numbers should be an analog of the Gaussian model 
with general D-brane insertions. Both are related by an open-closed duality in the spirit 
of ^ES]. 

Double-scaling limits of topological string theory have been studied before, mostly 
in the context of the Dijkgraaf-Vafa correspondence (see [26] for a recent example with 
references to the literature j^. In that case one studies type B topological string theory 

^[27] considers the double-scaling limit of Chern-Simons theory on the sphere, and recovers the c — 1 
behavior at the conifold. 
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on Calabi-Yau backgrounds related to the deformed conifold. The topological string 
amplitudes are given by matrix integrals, which can then be analyzed with the standard 
technology developed in the study of two-dimensional gravity. In contrast, we look at 
type A topological strings in backgrounds related to the resolved conifold, where the 
amplitudes are given by sums over partitions, and the connection to 2d gravity is more 
surprising (although it could have been suspected with some hindsight, due to the results 
of [an [221 [23]). 

The organization of this paper is as follows. In section 2 we review some results 
on phase transitions and phases of topological string theory. In section 3 we introduce 
perturbative topological string theory on Xp, write down its partition function, and relate 
it to various close cousins, in particular, we show that the model can be regarded in a 
precise sense as a q-deformation of Hurwitz theory. In section 4 we solve the model at 
genus zero (i.e. the planar theory) by using standard matrix model techniques. In section 
5 we extract the critical behavior of the planar theory and we show that including higher 
Casimirs leads to multicritical points a la Kazakov (albeit we do the detailed analysis 
only in the undeformed case, i.e. in Hurwitz theory). In section 6 we analyze the higher 
genus theory. To do this we make an ansatz for Fg{t) which generalizes the results of 
Goulden, Jackson and Vakil for Hurwitz theory [221 12H] to the deformed case. We also 
define the double-scaled theory and give strong evidence that it is described by 2d gravity. 
In section 7 we compare the structure found in the perturbative setting with the structure 
that one finds in the non-perturbative definition given by q-deformed 2d YM, and we 
speculate on the possible consequences for 2d gravity. In section 8 we point out evidence 
for a new open/closed duality generalizing the results of Aganagic and Vafa in [29]. In 
Appendix A we recall the very useful Lagrange inversion formula and apply it to various 
examples of the paper. In Appendix B we give some useful integrals for the saddle-point 
analysis of section 4. 

2 Phase transitions in topological string theory 

For simplicity, we will assume that the Calabi-Yau X has a single Kahler parameter t, 
i.e. h^'^{X) = 1. This is usually taken as a complexified parameter of the form 

t = r + i9, (2.1) 

but since we will be interested in questions regarding the convergence of the prepotential 
as a function of e~*, the 6 parameter is not relevant and we will set it to zero. 

2.1 Phase transitions in the planar free energy 

As we mentioned in the introduction, when t is large (in the so-called large radius regime) 
the geometry probed by string theory can be regarded as a classical geometry together 
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with stringy corrections. This is well reflected in the structure of the prepotential -Fo(^) 
or genus zero topological string amplitude, which in the large radius regime is of the form 



„ oo 

F,{t) = -t' + Y,No,ke-''. (2.2) 

k=l 

In this equation, C is the classical intersection number for the two-cycle whose size is 
measured by t. The infinite sum in the r.h.s. is given by worldsheet instanton corrections, 
which are obtained by "counting" (in an appropriate sense) holomorphic maps from 
to X. The instanton counting numbers A^o.fe are genus zero Gromov-Witten invariants, 
and we have chosen units in which ig = v^27r. 

The series of worldsheet instanton corrections, regarded as a power series in e~*, has in 
general a finite radius of convergence which can be obtained by looking at the asymptotic 
growth with k of the numbers A^o.fc- We will characterize this asymptotic growth by tc 
and by a critical exponent 7: 

No,k ~ k^oo. (2.3) 

When this holds, the prepotential behaves near tc as 

Fo(t)~(e-*^-e-*)2-^ (2.4) 

It turns out that typical Gromov-Witten invariants of Calabi-Yau manifolds behave 
asymptotically as 

No,k ~ Tjf^, k^oo. (2.5) 
k'^ log k 

This is of the form fl2.3p . with critical exponent 

7 = (2.6) 

and subleading log corrections. This behavior was first established in PU] in the example 
of the quintic, and since then it has been verified in other examples, like for example in 
local P^, where the critical radius is given by 



= —-^--ReGf -,-,!;!) - 2.90759 (2.7) 



r,3 



3'3' 



and G is the Meijer function. 

The subleading log in (12.51) leads to log corrections near the critical point (also referred 
to as scaling violations) of the form 

Fo{t) ~ (e-*= - e-*)'log(e-*= - e"*). (2.8) 
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This is the genus zero free energy of the c = 1 string at the self-dual radius, once the 
scaling variable e~*'= — e~* is identified with the cosmological constant [3],|1]. 

The behavior of the prepotential gives a precise quantitative meaning to the distinction 
between classical and quantum geometry. We will refer to the divergence of the large 
radius expansion at t = as a phase transition with a critical exponent 7 defined in 
( I2.5p . The phase with 

t > tc (2.9) 

where the expansion (12.21) is convergent, is called the large radius or Calabi-Yau phase, 
where classical geometry makes sense (albeit it is corrected by worldsheet instantons). 
When t < tc, the nonlinear sigma model approach is not well defined, and classical 
geometric intuition is misleading. 

There are two possible approaches to this problem. In [311 [33] , Aspinwall, Greene and 
Morrison proposed to study this small area phase by using mirror symmetry. One finds 
that the point t = tc separating the two phases corresponds to the discriminant locus of 
the mirror moduli space (the conifold point). For t < tc, mirror symmetry makes possible 
to compute the prepotential by analytic continuation. The second approach is due to 
Witten, who proposed in [34j the linear a model description of A/" = 2 sigma models 
on Calabi-Yau manifolds largely as a way to understand this phase transition. In this 
formalism, one recovers the mirror symmetric description, and in addition the stringy or 
non-geometric phase can be described in a precise way. It is typically a Landau-Ginzburg 
orbifold perturbed by twist operators. 

2.2 Higher genus 

In order to describe the phase structure of the model we have relied on the behavior of 
the prepotential, i.e. the planar free energy. It is natural to ask what happens when 
higher genus topological string amplitudes are taken into account. These amplitudes, 
which we denote by Fg{t), can be expressed in the large radius limit in terms of genus g 
Gromov- Witten invariants Ng^k'- 

00 

Fait) = J2^9,ke-'''. (2.10) 

k=l 

Here we have omitted contributions from degree zero maps as well as classical pieces for 
g = 1, which will not be important in the discussion. It turns out that the Gromov-Witten 
invariants have the asymptotic behavior [3] 

Ng,k ~ A;(^-2)(i-g)-igH,^ (2.11) 

where t^ is the critical radius obtained at genus zero and it is common to all g, and 7 is 
the critical exponent that appears in fl2.3l) . This is equivalent to the following behavior 
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near the critical point 

~ci log (e-*^-e-*; 



F,(t)~c,(e-*=-e-*)(i-^)(2-7)^ g>2. 



(2.12) 



In conventional topological string theory, as we have mentioned, 7 = 0, but the more 
general form we have written above will be useful later. 

We then see that the phase transition ait = tc is common for all Fg{t), and the critical 
exponent changes with the genus in the way prescribed by (12. lip . This sort of coherent 
behavior in the genus expansion is not obvious, but seems to characterize a wide variety 
of systems that admit a genus expansion (like for example matrix models, see [35] for a 
review). When this is the case, one can define a double-scaling limit [Hllinilll] as follows. 
Let us consider the total free energy F as a perturbative expansion in powers of the string 
coupling constant Qs'. 

00 

F{g,,t) = J2m9'r'. (2.13) 

9=0 

We define the double-scaled string coupling as 

K = ag,{e-'^-e-Y'-\ (2.14) 

where a is an appropriate constant. We can then consider the limit 

t tc, (?s — ^ 0, K fixed. (2.15) 

In this limit, only the most singular part of Fg{t) survives at each genus, and the total 
free energy becomes the double-scaled free energy 

FM = fo^~' + /i log + 5^ fgK'^~\ (2.16) 

where fg = a^~'^^Cg. It is also customary to express the double-scaled free energy in terms 
of the scaling variable z = 

It turns out that, in some cases, one can determine the coefficients fg in closed form. In 
the double-scaling limit of matrix models, they are governed by a differential equation of 
the Painleve type [35j . In the case of topological string theory on Calabi-Yau manifolds, 
it was conjectured in |4j that, in terms of a natural coordinate 

/i~e"*^-e-* (2.17) 

which in the mirror model measures the distance to the conifold point n = 0, the double- 
scaled free energy is universal and reads 

1 1 °° R 

FM = log ^ - - log . + g ^^i^/-'-'- 

This is exactly the all genus free energy of the c = 1 string at the self-dual radius (for a 
review, see [36]). This behavior has been checked in many examples [32| [37]. 
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3 Topological string theory on bundles over 
3.1 The model and its partition function 

In this paper we study perturbative topological string theory on a Calabi-Yau threefold 
given by the total space of a rank two holomorphic complex bundle over P^. It is a well- 
known result that these bundles split into a direct sum of line bundles, therefore these 
spaces are of the form 

0{p - 2) ® 0{-p) ^ F\ (3.1) 
where p G Z. There is an obvious symmetry 

+ 2 (3.2) 

and we can restrict to p > 0. 

The topological string partition function on the space Xp is the exponential of the 
total free energy (I2.13P 

=exp Fxp{gs,t) (3.3) 

where 

oo oo oo 

9=0 9=0 d=l 

and Nd,g{p) are the Gromov-Witten invariants. This partition function can be com- 
puted by using the topological vertex [38] , and the same result is obtained from the local 
Gromov-Witten theory of curves of |S] (see section 2 of [6j for a useful summary). In order 
to write down the result, we recall some necessary ingredients. First of all, we define the 
q-number [n] as 

[n] = g"/2 - /2, q = e-9\ (3.5) 

A representation RoiU {oo) can be represented by a Young tableau, labeled by the lengths 
of its rows {li}. The quantity 

^{R) = Y,h (3-6) 

i 

is the total number of boxes of the tableau. Notice that a tableau R with i{R) boxes can 
be regarded as a representation of the symmetric group Su^r) of i'(-R) elements. Another 
important quantity associated to a tableau is 

«:H = ^/i(/,-2i + l), (3.7) 

i 

which is closely related to the second Casimir of R regarded as a U{N) representation. 
Finally, we introduce the quantity 

Wn = q-^^" n ru L ,1 - (3-8) 
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This is in fact a specialization of the topological vertex. The product in the r.h.s. is over 
all the boxes in the tableau, and hook(n) denotes the hook length of a given box in the 
tableau. 

The topological string partition function on Xp is given by 

= J2 WnWR^q^P-'^^^/^-lf^'^Pe-'^''^' (3.9) 

R 

where i?* denotes the transpose tableau (i.e. the tableau where we exchange rows and 
columns). The above expression has to be understood, as in [38], as a power series in e~*. 
At every degree it gives the all-genus answer. One can easily compute the first few terms 
in the expansion: 

F^^it) = i-iye-' + ^(2p2 - 4p + l)e-2* + -^^(1 - 6p + 3p2)(2 -6p + 3p'')e-^' 

8 54 



+ C(e 



-4t\ 



(-l)p 

+ ^^^(-2 + 14p - 19/ - 20 + 45 - 24 + 4p'')e-^* + 0(e-^*), 

(3.10) 

and so on. The Gromov-Witten invariants Ng^d{p) are polynomials in p of degree 2d—2+2g 
with rational coefficients: 

2d-2+2g 

N,Ap)= E ^l^P'- (3-11) 

As we will see, the coefficient of the highest power in (13. lip has a simple geometric 
interpretation in terms of classical Hurwitz theory. Note that, due to the definition we 
are using of Qs, the Gromov-Witten invariants N^^g differ in a sign (—1)^'^ from the 
standard ones. This is of no consequence for most of our discussion, but the extra sign 
should be taken into account when extracting integer Gopakumar-Vafa invariants Ug^p) 
[39] from the above expressions. 

The partition function (13. 9p on Xp is calculated, as explained in [Sj, by considering 
the equivariant Gromov-Witten theory on Xp with respect to scalings of the line bundles 
over P^ and by using the antidiagonal action. This is, therefore, an equivariant partition 
function, and leads to a topological string theory which is different from the standard 
ones. For p = 1 this theory specializes to (standard) topological string theory on the 
resolved conifold. For p = 2, the geometry is that of C x Ai. The standard topological 
string theory has a trivial partition function on this space, while the equivariant theory 
considered here has = Zx^. See [ID] for more details on the equivariant theory from 
the point of view of mirror symmetry. 

Another way to understand the geometric content of this non-standard theory is to 
consider the equivariant version of topological strings presented in section 7.6 of [H]. Let 
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us look at equivariant topological string theory on local P^. This model contains, on top 
of the Kahler parameter t, three equivariant parameters Xj, i = 1,2, 3, related to twisted 
masses in the sigma model. The partition function depends now on the Xi, 

Zp2{xi,X2,X3), (3.12) 

and can be computed again by using the topological vertex (explicit formulae are given 
in [S]). It is completely symmetric under permutations of the and, if we set Xi = 1 
we recover the standard local partition function. One can deduce from the explicit 
formulae in [H] that 

Zp2(l,0,0) = Z^3. (3.13) 

This is easy to understand, since the hyperplane class of P^ in the total geometry is indeed 
a P^ with normal bundle (9(— 3)©(9(1). By taking X2 = X3 = we are considering the local 
geometry of this curve inside C(— 3) P^. For p > 3 we can regard Zx^ as a specialization 
of equivariant topological string theory on more complicated toric backgrounds. 

Finally, we point out that (13.91) can be interpreted in terms of representations of Hecke 
algebras by using the results of [12]. 



3.2 Relation to Hurwitz theory 

Before analyzing this partition sum in detail, let us relate it to other generating func- 
tionals. The quantity Wr turns out to be a q-deformation of the dimension dn of the 
representation R of S'f(K). As gs — > 0, one has that 



d 



\i{R)\\ ^ ^ 

This suggests taking the following limit, 

(7s ^ 0, t ^ 00, p ^ 00, (3.15) 

in such a way that 

P9s = r^/N, (-l)Pe-* = {gsNfe--\ (3.16) 

and Ti, T2 and N are new parameters that are kept fixed. It will be convenient in what 
follows to introduce a 't Hooft parameter as 

T = gsN. (3.17) 

In the limit (I3.15p - (l3.16p one has T ^ 0. The partition function becomes: 

Zx, - Zhu™.. = E(^)'^''^''^^'"^"^'''^""'^''^- (3-18) 

R 
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This is the generating functional of simple Hurwitz numbers of at all genus and degrees. 
In order to see this, it is useful to recall some basic ingredients of Hurwitz theory. 

Hurwitz theory studies branched covers of Riemann surfaces. The structure of the 
covering will be labeled by a vector k of nonnegative entries, and we define 

= 1^1 = J] z^.^Hf^kjl (3.19) 

j 3 3 

If a cover has degree d, one necessarily has 

lik) = d. (3.20) 
Given m branched points on P^, their branching structure can be specified by m vectors 

(3.21) 

The Riemann-Hurwitz formula relates the entries of these vectors, the degree of the cover 
and the genus g of the covering surface as follows: 

2g-2 + 2d = ^{t{k^'^) - \k^'^\). (3.22) 

i 

The number of Hurwitz covers at degree d is given by the classical formula 

(, s 2 m 

In this equation, 

Uk) ^ &1M(£M, (3.24) 
dR 

— * — * 

where C{k) is the conjugacy class associated to the vector k in the symmetric group of 

— * — * 

£{k) elements, and XR{C{k)) is its character. 

A branched point whose branching vector is of the form 

ki = id-2,l,---) = {2) (3.25) 

is called a simple branched point. A particularly important example of Hurwitz theory 
occurs when we have one non-simple branched point specified by a vector k, and m{g, k) 
simple branched points. By Riemann-Hurwitz, the number of simple branched points is 
given by 

m{g,k) ^2g-2 + \k\+e{k), (3.26) 

— * 

where we used that d — £{k). Since 

fR{d- 2,1,0, ■■■) = Ikr, (3.27) 
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The Hurwitz number for this configuration is given by 

<i(^)= E (3-28) 

£(R)=d^ ^ ''' 

The case where there are only m simple branched points is a particular case of the 
above. It is obtained when the branching structure is trivial, 

k = [d,^,---) = \^, (3.29) 

and 

m = 2^ - 2 + 2rf. (3.30) 

The resulting Hurwitz number is called a simple Hurwitz number, and it is given at genus 
g and degree d by 

E (^^I'^'^^Z^)''"''''' (3-31) 

where the sum is over representations R with fixed number of boxes equal to the degree 
d. 

We can now rewrite fl3.18p as 



2,-2+2. (3-32) 
9,fl!V ^ 2o-2 + 2d ' 

where in the second line we have traded the sum over m by a sum over and we used 
(13.301) . The model described by 03.181) has been studied in detail due to its connection to 
Hurwitz theory. From the physical point of view, it was analyzed in pil HSl Hlj, and in 
the mathematical literature it has been studied for example in [22l [28] . 

The free energy of ^Hurwitz describes connected, simple Hurwitz numbers Hj^^{l'^)': 

2g-2+2d 

FHu.w.tz = l0g^Hurw.t. = ^ ^'"'^ E ^'^''K^i^'T ^Oo- 2 + 2dV ^^'^^^ 

g>0 d>0 I y + J- 

If we compare this to the total free energy written in 03.41) in terms of Gromov-Witten 
invariants, and take the limit 03.15p - 03.16p . we find 

lim p-^»--(-l)^iV,,(p) = pg^fl;,), (3.34) 
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The l.h.s. is precisely the coefficient of the highest power in p of dmia. We can therefore 
interpret the Gromov-Witten invariants of this model as q-deformed connected, simple 
Hurwitz numbers, since they promote to polynomials of degree 2g — 2 + 2d 

(which is equal to the number of simple branch points). 

We also point out that a closely related model to the partition function Zx^ and 
its Hurwitz limit is four-dimensional M = 2 supersymmetric gauge theory with 

Casimir operators, considered in [15] from the point of view of Nekrasov's instanton 
counting [l6]. When only the ffist and the second Casimir operators are turned on, this 
model is equivalent to the Hurwitz model fl3.18l) . and via the Hurwitz/ Gromov-Witten 
correspondence of ^S], it is related to topological string theory on P^. The q-deformed 
model Zxp we are considering here is in turn related to the five-dimensional or K-theoretic 
version of the gauge theory (see for example [H]). 



3.3 Relation to open Gromov-Witten invariants 

There is a generating functional in open Gromov-Witten theory which is also closely re- 
lated to Hurwitz numbers. Consider a Lagrangian submanifold with the topology C x 
in C^, and with framing /. The generating functional of open Gromov-Witten invari- 
ants can be computed with the Chern-Simons invariants of the unknot [50] or with the 
topological vertex [3SIEI], and reads: 

Z{V) = iy,jg^"«/2(-l)'(^)^Tr^ V, (3.35) 

R 

where is a [/(oo) matrix source of open string moduli that takes into account the open 
string sectors associated to different winding numbers. Let us now rescale V — > e~*K and 
take the hmit (appropriate for the open sector) 

Qs^^, t ^ oo, / ^ oo, (3.36) 

keeping fixed 

fgs = T-s/iV, (-l)^e-* = gsNe--\ (3.37) 
The generating functional fl3.35p becomes 



R 

\2g-2 + \k\+m)\ 



(3.38) 



^This result, for the genus zero case, was already derived in [47] 
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where 

T^{V) = Y[{TiV^Y\ (3.39) 

i 

f l3.38p is a generating functional for the more general Hurwitz numbers (13.281) . If we now 
consider the free energy associated to the open Gromov-Witten functional, and we write 
them in terms of open Gromov-Witten invariants F^j:{f), 

oo 
k 9=0 

we find that 

f^°° {2g - 2 + \k\ + £{k))\ 



This expresses the more general Hurwitz numbers (I3.28P as limits of open Gromov-Witten 
invariants. We point out that the open Gromov-Witten invariants ^(/) can be ex- 
pressed in terms of triple Hodge integrals [52| [53] , while the Hurwitz numbers (13.281) can 
be written in terms of simple Hodge integrals [51] , and the relation (13.411) has been noted 
in [551 EE] by considering their integral expression. 

We now relate the closed Gromov-Witten invariants Ng^d{p) in the background Xp to 
the open Gromov-Witten invariants F^^{f). It is clear that (I3.35P reduces to (13.90 if we 
set / = p — 1 and if we give a value to the source V in such a way that 

Tr^ V = {-if^^^e-^^^^'WRt. (3.42) 

There is indeed a choice of the (infinite-dimensional) matrix V which produces this, with 
eigenvalues Xi = g*~^/^e~*. With this choice, one has that 

rdV) = i-lfe-'^'^' n , . ^ (3-43) 

• [q2 — q 2)i^j 



If we expand this in powers of g. 



Si 



oo 



T,K\/)=e-^(^)*5^^;l'=l+X,^, (3.44) 

3=0 



where g are rational numbers, we find 



NM= E (3.45) 

e{k)=d 9'=o 
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which can be used to express the q-deformed Hurwitz number Ng^d{p) in terms of triple 
Hodge integrals. 

There is a more elegant way to relate the closed and the open invariants, by using the 
underlying integer invariants. It was shown in section 5 of [?T] that, when closed string 
functionals are obtained from open string functionals by "condensing" the source term as 
in (13.421) . the Gopakumar-Vafa invariants ng^d of the closed geometry can be written in 
terms of the open BPS invariants introduced in [50l [57] . In the case of framed branes in 
C^, the integer invariants Nii^g{f) depend only on the genus g and a representation R of 
f/(oo). It is easy to obtain from the results in [?T] that 

UgAp) = N^_^,gip - 1). (3.46) 

d boxes 

The tableau in the r.h.s. corresponds to the trivial representation of Sd, and in deriving 
this relation we have used the conventions for open BPS invariants in [531 EZj- 



4 The planar limit 

The explicit expression (13. 9p gives the total partition function but it is not useful in 
extracting the functions Fg ^(t) in a compact way. If we are interested for example in the 
asymptotic behavior of the Gromov-Witten invariants for large degree, we would like to 
have closed expressions for the free energies at fixed genus but at all degrees. These are 
the kind of expressions that can be obtained by using mirror symmetry. Unfortunately, 
there is no complete mirror description of the Xp backgrounds that can be used efficiently 
to compute Fg ^(t). Some results concerning this description have been recently obtained 
in [lO] , but most of the general ingredients (like Picard-Fuchs equations) are still lacking. 

In these circumstances, we have to extract the Fg ^{t) directly from the sum over 
partitions (13.91) . In this section we solve the planar limit, i.e. we find the genus zero 
contribution or prepotential, by doing a saddle-point analysis of (13. 9p . Some aspects of 
our analysis were worked out in [Ti] . 

4.1 Saddle— point analysis 

The sum over partitions that defines Zxp in (13. 9p can be rewritten as 

= Y,w'^ g(p-2)««/2 ^-mt_ (4.1) 

R 

where for simplicity we have absorbed the sign (— 1)^(^)p in the Kahler parameter t. To 
analyze this in the saddle-point limit, we will adopt a strategy used by Kostov, Staudacher 
and Wynter in |2lj to analyze the partition function of Hurwitz theory (13.180 . We will 
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introduce a "fake" dependence in the theory, which allows for standard large N analysis, 
and we will then extract the genus zero result. 

To see this, we first notice that (14.11) admits an evocative representation in terms of 
a q-deformed group theoretical quantity of U{N). Let {Zj} be the lengths of rows in a 
Young tableau, and let hi = li + N — i. Consider the q-deformed quantity 

and the familiar quantum dimension of an irreducible representation R 

l<i<j<N J 

The brackets denote as usual q-numbers as in (13. 5p . Using the equality 

dim,i?^^_,^/,^^^ (4.4) 



q^R 

we can write 

Zx, = J2f^^yq(p-^)^n/2^-^m_ (4.5) 

\ q R / 

However, the actual equality in (14. 5 p holds only when the r.h.s is expanded as an asymp- 
totic series in this immediately suggests to exploit the general large N techniques 
in investigating (14. 5p . 

The planar theory can be in fact immediately analyzed along the usual strategy orig- 
inally introduced by Douglas and Kazakov for QCD2 [58j. We introduce the auxiliary 't 
Hooft parameter T = Nqs as in (13.171) and continuous variables in the standard way: 

I = I - ^ + 1 - ^(x) -x + l = h{x), (4.6) 

The delicate point is to evaluate the large limit of the deformed measure. The numer- 
ator of qflR leads to 

N hi N hi , . 

lognn(^'^ - q-^^f = 2EEl°g 2sinh^.^, (4.7) 

1=1 j=l i=l j=l 

which becomes in the large limit 

1^ dx dy\og2smhT ^ =^ dx\^-^ - - + U,{e-'''^)j. (4.8) 
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The denominator of qQpi cancels against the denominator of (dimgi?)^. The numerator 
of the quantum dimension leads to a sinh (x — y) interaction, as explained in [121 [ISl [H] 
in a related context. Then we can write the effective action controlling the leading large 
contribution as follows 



1 .1 2 



S = — / dxdylog 
'0 Jo 



2 sinh —(h(x) — h{y)) 



T 



2 . (4-9) 



+ ^ dxh{x){t -{p- 1)T) + ^ ^ dxh'ix) + {p-l)l-^- It. 

The planar theory can be, thus, understood as coming from a matrix model: we notice, 
in fact, that the effective action can be derived from a Chern-Simons-like matrix model 
[59] with a potential V{h) of the form 

Vih) = |Li2(e-^^) + (t - (p - mh + ^/^^ (4.10) 
and the saddle-point equation is simply 

I dh'pih') coth -h')=ph + ^ log(l - e-™) + 1 - (p - 1), (4.11) 
where the density p{h) is defined in terms of the inverse function x{h) as follows 

(4.1, 

Because of the posit ivity constraint 

hi> h2> ■■■hN>0^ h{x) >0, (4.13) 

which the Young tableaux variables hi must satisfy, the support of p{h) will be chosen in 
the interval [0, a]. 

The above equation can be related to a standard Riemann-Hilbert problem: to this 
aim we first introduce x = 1 — hT, z = 1 — h'T, and then we pass to exponential variables 

s = e^ y = e\ (4.14) 

In terms of these variables the saddle-point equation reads 

plogs-(t + l) + (p-l)(T-l)-21og(l-e-^s), (4.15) 



Je-i3 y 



s-y 

with —(3 = 1 — Ta. The normalization of p is now 



dy(^ = T. (4.16) 

y 
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The support of the density p{s) comes from the original tableau variables h, i.e. [0, a]. 

To solve the saddle point equation fl4.15p we need a further ingredient, namely we have 
to choose an ansatz for the density p{s). In order to recover the large radius expansion 
in e~* the analogy with QCD2 suggests to choose a chiral, one-cut ansatz: for x G 
[— /3, —7] (—7 < 1) the p(s) is arbitrary, while for x e [—7, 1] we require p to be equal to 
1. With this assumption, we easily arrive at our final form for the saddle-point equation 

= — logs log(l-e^s). (4.17) 

y s — y 2s 2s s 

Notice that in principle this equation and its solution depend on the parameter T = 
Nqs, therefore the genus zero free energy of this model Fo(T, p, t) depends on T, as well 
as on p and t. However, the parameter N was introduced by hand and does not appear 
in the original model. Consistency of our procedure requires then 

Fo(r,p,t) = ^FoiT=l,p,t) = ^X'ip.t). (4.18) 

The reason for this is that N^Fo{T,p,t) should equal g;^FQ^{p,t). We can indeed verify 
this in detail as follows. If we perform the changes of variable 

s i-^ z = se^'^ and y ^ v = ye^~^, (4.19) 

accompanied with the redefinitions 

7 = 7 + 1- T, P = P + 1-T and p{v) = p{ve^-^), (4.20) 

the saddle point equation takes the form 

"-^ = f log.-f-ilogd-e''.). (4.21) 

V z — V 2z 2z z 

The normalization condition instead reduces to 



-7 



—p{v) = -7. (4.22) 

We can now verify fl4.18p by computing the T— dependence of the derivative of Fq with 
respect to t. In terms of the original Young tableaux variables, we have 

OF, 



dt 



[ p{h)hdh--. (4.23) 
Jo 2 



Now if we set /i = 1 — log(fe^ with the help of the redefinitions (14.191) and fl4.20p . 

the above equation can be cast in the following form 



18 



as needed. In the following we shall set T = 1, since in this way the planar free energy of 
the above matrix model equals the prepotential Fq ^{p, t); the T dependence is eventually 
recovered through the relations fl4.19p . (14.201) and (14.241) . 



4.2 Solving the saddle— point equation 

The solution to the integral equation (I4.17P can be written in terms of the effective resol- 
vent function 

Je-/3 V Z-V 

which is given by Muskhelishvili-Migdal's formula 

dv I log - I - log (1 - ve^ 



27rz^^ 'Jc{z-v)v ^{v-e~P){v-e~'<) ^ ' 

where the closed contour C encircles the support [e~^,e~'>'] of the distribution p{v) with 
counterclockwise orientation in the complex z-plane. If we choose the square root and 
logarithmic^ branch cuts in (I4.26P as indicated in Fig. [H then since the integrand decays 




Figure 1: The contour Ci surrounds the branch cut (— oo,0] of log(z), C2 encircles the 
cut , 00) of log (1 — e'^'u), while C encloses the physical branch cut [e~^,e~'>'] . 

as at |f I 00, we can deform the contour of integration C so that it encircles the 
cuts of the two logarithms. This deformation picks up an additional contribution from 

^We have always defined the logarithm function as having a branch cut along the negative axis. This 
choice implies that log(l — e^w) has a cut for w > e^^ . 
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the pole at w = z and we find 



V 



e ■2\/z — e 



-y + 6^2/ 



z — e 



\Jz - e-^ + ^/z- e-^ 



z 



( ( ^ y^7^ V\ 

(e-7 - e-/5) 



£+7 
e 2 



v/(;z-e"T)(z-e-/5) plog 



62+62 



log 



_2 , 

6 2+62 

6 2 — e 2 



1_ 1 Le(/5+7)/2 

2 2z 2z 



^(^_e-7)(^-e-/3). 



(4.27) 



Because of the chiral ansatz, the boundary condition for large z explicitly depends on the 
extremes of the interval and it imposes that 1-1^(2;) ~ —l/z. Then the vanishing of the 
constant term in (14.271) at infinity implies that 



[l-p) log 



i. 1 

62+62 



log 



62—62 



t 

2' 



(4.28) 



The matching of the subleading term in the expansion {z ^) instead imposes the following 
endpoint equation 



(p-l) logl 



6 2+62 



log 



6 2—62 



(4.29) 



Our goal, in the following, is to solve these two equations in a closed form. The first step 
is to show that (14.281) and (14.291) can be reduced to a polynomial equation for an auxiliary 
unknown w. First, we introduce two intermediate variables x and y defined as 



X 



6 2+62 



_2 -i. 
6 2—62 



(4.30) 



which, because of their definition, have to obey the following inequalities: a; > 0, y > 
0, X — y = 6~2 > 0. The saddle point equations (I4.28P and (14.291) now read 

-t/2 l-p ( 2 _ 2\P t/2 p-l 

f ^ ^ IJ_ _ 1 = 1 = 0. (4.31) 

y y 

The second equation can be easily solved with respect to y {y = e~*/'^x^~^) and, after a 
trivial algebraic manipulation, we are left with 



X 



^ - 1) = 6- 
20 



(4.32) 



The form of (14.321) suggests that the natural variable for our problem is 

w = X p , (4.33) 
in terms of which fl4.32p becomes a polynomial equation 

e-* = w^P~^^"-' - w^P-'^' = f{w) . (4.34) 

Because of the inequalities constraining x and y, we are interested in the solutions of 

_ t 

(14.341) that satisfy w > e . Once we have solved (I4.34p . the endpoints are then 
recovered through the relation 
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-2 log (w-P/^ + e^*/2w-5(P-i)pj and (3 = -2 log (w~p/^ - e"*/2u7-5(p-i)p j . (4.35) 



The existence and the properties of such solutions can be investigated by studying the 
behavior of the r.h.s. of fl4.34p and, in particular, it is sufficient to focus on the region 
< w < 1, since e~* belongs to the interval (0, 1] when t > 0. In this region, for p > 2, 
the derivative of the r.h.s. of (14.341) . i.e. 

f{w) = w^P'^^'~\{p - 2)p - (p - Ifw), (4.36) 

has a unique zero given by 

PjP - 2) ^ ^ , . 

= (p _ 1)2 < 1- (4-37) 

Moreover it is positive for < w < Wc and negative for w > Wc. 

Then solutions of fl4.34p exist if and only if e~* is less than the maximum located at 
w = Wc. In other words, we must have e~* < f{wc). This inequality is saturated by a 
minimum value of t given by 

tc = log {{p{p - 2)y^'~P\p - l)2(p-i)^) . (4.38) 

Above this critical value, there are two solutions: the former is smoothly connected to 

w = 1 while the latter is smoothly connected to w = 0, as fig. [2] shows. 

t 

Because of the bound w > e p'-p-^) , the solution which appears to describe the regime for 
large t is the one near w = 1. This solution can be constructed as a series in e~*: 



w = 1+ y x„,e 



n=l 



The coefficients Xn are determined, in a closed form, by means of the Lagrange inversion 
formula applied to fl4.34p (see appendix A). We obtain, after reinserting the sign factor 
through e~* — >■ (— l)^e~*. 



oo n—1 

1 + 5]; n + ^ - - 1)') — (4-39) 



w 



n=l k=l 
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0.2 



0.4 



1 



Figure 2: Plot of the r.h.s. of (14.341) in the region < w < 1. The straight hne represents 
the constant value e~*. 



An intriguing feature of this expansion, even if it is not evident, is that its coefficients are 
integer numbers for integer p. This suggests that the relation between w and e~* is a mirror 
map, which is indeed characterized (for yet unknwon reasons) by integer coefficients. 
It is worth noticing that the above series can be also summed. Our equation belongs, in 
fact, to a special class of algebraic equations known as trinomial equations. These are 
equations of the form 

- aw" + 6 = 0, with n > s, (4.40) 

and it was shown by Birkeland (see Appendix A in [60]) that they can be solved in terms 
of generalized hypergeometric functions. In our case, we find 



w 



1 n — 
- + 

n n 



1 

n-1 



c 



(4.41) 



-{t-tc) 



with n = {p — 1)^ and C = ^ 
It is instructive to investigate the radius of convergence of the series (14.391) . The simplest 
way of finding it is to look directly at (I4.4ip . From the theory of hypergeometric series, 
we know that our series converges only if C ^ namely t > t^- it starts diverging for the 
same value t^ for which the solutions to the endpoints equations stop existing. The same 
result can be reached by an estimate of the asymptotic behavior of the Xn through the 
Stirling formula. As we will see in the next section, this is also the radius of convergence 
of the prepotential. 

In this analysis, we have neglected the cases p = 1,2, which deserve special attention. 
In these two cases, solutions exist for any value of the parameter t greater than zero and 
the solution can be given in terms of elementary functions: fl4.34p is, in fact, of the first 
order in w both for p = 1 and p = 2. 
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At this point we should remark that the above solution can be potentially shadowed 
by the breaking of the chiral ansatz. In fact, the endpoint 7 may reach —1 {h = 0) for 
a value tc > t^. However, there are various reasons why this transition is not relevant 
for our analysis. First of all, this transition depends on the value of T: it exists only for 
small value of T and when T > 1.68, the transition disappears completely. This indicates 
that the transition is not intrinsic to our original model (14.11) . but rather an artifact of the 
matrix model description. Second, this transition is related to the Gross- Witten-Wadia 
transition flbi [T6] , therefore it does not characterize the expansion (14. 5 p in itself. The 
transition appears only when this expansion is regarded as the strong coupling expansion 
of a unitary matrix model. 



4.3 The planar free energy 

When the 7) are chosen so that equations (I4.28P and (I4.29P are satisfied, the form of 
the resolvent greatly simplifies. In fact, only the logarithmic parts survive and we are left 
with 



uj{z) 



z 2z 



+ 



P 1 



e 2 v^^^^"e^ + 62 y/ z — e^f^^ 



( 



\ 



) 



+ -lo^ 

z 



(e-7 - e-/3) 



v 



(4.42) 



^/z - e-T + ^z - e-^ ' 

Its jump around the cut (e~^,e~'^) determines our genus zero density function p. We get 
p(s) = arctan | \\ — ) H — arctan | \ \ — ^ | — — arctan | e2~2 



TT 



e-T - s 
s — e^^ 



J 



p e ^ — s 



e — s j 71 \ Vs — 

(4.43) 

We now possess all the necessary ingredients for computing the partition function at genus 
zero. We proceed as usual and we first evaluate its derivative with respect to t. It is given 
by 



d{Fo{t,p)) 
dt 



' ds 7 
— logsp(s) - — 
-0 s z 



1 r ds {log sf{p{s)y, (4.44) 

2 Je-fl 



where the contribution —7^/2 comes from the region where the original density is constant. 
The last integral in (I4.44p can be computed in a closed form and only in terms of w. To 
this purpose, we shall first change the integration variable from s to 



(e ^ — e ^) cos( 



(4.45) 
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The integral then reduces to 



d0 



2-2 

in 



_/3+7 
e ^ p 



27r( (e-Z'-e-T) cos (/)+e-''+e-T) 



-'3+e-T-(e-T-e-'3) < 



(4.46) 



and it can be exactly performed as shown in appendix B and the result in terms of the 
two endpoints is 



I 



-Li2 



6/3/2 + e7/2 



+ (p - 2) log2 



-/3/2 



+ e- 



-7/2 



— p log^ 



;/3/2 _^ g7/2 



By employing the expression of P and 7 in terms of w (I4.35p . we obtain 



^ = ^P^P - 2) log^(u^) - Li2(l - w) 



(4.47) 



It is interesting to notice that taking another derivative of the partition function with 
respect to t produces an even simpler result. Since (14.341) fixes the derivative with respect 
to t of If to be 

// N (W— 1)W / , , X 

^ W = -7 TT. 7 (4-48) 



(p — l)'^w — {p — 2)p' 



an easy computation shows that 



log(w). 



(4.49) 



To obtain the actual free energy of genus zero, we have to integrate (14.471) over t. This 
integration, however, can be transformed in an integration over w with help of (14.481) . We 
can write 



dw 



nj,p) 



p{p-2) ^ 1 



w 



w — 1 



p{p - 2) \o^{w) - Li2(l \ , 



(4.50) 



which, integrated over w, gives 
Fo(w(t),p) = \oi{w) + Li3(l - w)^p{j> - 2) log(w)Li2(l - w) 



+ 



6 

{P - 2)p 



(log(l - w) \oi[w) + 2Li2(M;) log(M;) - 2Li3(w) + 2C(3)) 



(4.51) 



up to p— dependent constant. With the help of the polylogarithmic identities that con- 
nect polylogarithms of different arguments (see appendix B), we can reduce the partition 
function to its final form 



1 



P , 



Fo{w{t),p) =p{p - 2)Li3 1 - - + (p - l)^Li3(l -w)-^{p-2){p- ly \og'{w). 



6 



(4.52) 
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We can also provide a closed expansion for the prepotential Fq as a series in e~*. This 
is better done by working out the expansion of log w through Lagrange inversion and 
integrating (14.491) twice. In this way we obtain 

It is now easy to check the consistency of this result with the direct expansion of free- 
energy in (13.101) . justifying therefore a posteriori the choice of the chiral ansatz. We notice 
that the partition function for p = 1,2 greatly simplifies and we obtain, in both cases, 
the partition function of the resolved conifold, as expected, and up to a global sign due 
to our choice of gs, 

Fo^^(t) = -Li3 (e-*) , Fo^^(t) = Li3(e-*). (4.54) 

Finally, we note that the exact expression (14.531) only depends on p through (p — 1)^, in 
accord with the symmetry (13.21) . 



4.4 Comparison with Hurwitz theory 

As we mentioned in section 3, the model (14. ip can be regarded as a q-deformation of 
Hurwitz theory, which should be recovered in the limit (13.161) . This limit can be taken 
order by order in the genus expansion, as one can easily see by writing it in terms of the 
't Hooft parameter T. We have to take 

T^O, t^+oo, p^oo, (4.55) 

in such a way that 

pT = T2, (-l)^'e-*t = TV'S (4.56) 
are kept fixed. If we write the free energy of Hurwitz theory (13.331) as 

2-23 



E - ^"(Z^)' (4-57) 



^ Hurwitz 

9=0 

where 

^^ = rle-^\ (4.58) 

then in the limit (I4.56P one has 

p2-2gpX, ^ ^Hurwitz _ gg^ 

The planar limit of the Hurwitz model was analyzed in [21] by using matrix model 
techniques. It is easy to see that all of their results can be recovered from the planar 
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solution of the deformed model. In the solution of [21], a crucial role is played by an 
auxiliary variable Xi which is related to the endpoints [6, a] of the Young tableaux density 
through 

X = ri + 21og^. (4.60) 
The endpoint equations of [21] lead to an equation relating x to the parameter 

/z = xe-^. (4.61) 
The solution x(/i) to this equation is provided by Lambert's W function 

oo , k—1 
k=l 

which has convergence radius Hc = e^^. It follows from the results of ^21] that the 
prepotential F^'^^'^^*'^ {fj,) of Hurwitz theory satisfies 

^o"(/^) = X, (4.63) 
and from here one can find the power series expansion for the free-energy 

oo , k—3 

Fo"^(^) = 5^— (4.64) 

k=l 

Comparing the undeformed case with the deformed one, we realize that the variable 
X is playing here the role of the variable w introduced in (14.331) . In fact in the limit 
fH35D - (lt56D one has 

w-1^-4- (4.65) 

Also, the equation (14.640 is the limit of the equation (14.490 . and using this fact or the 
explicit expansions ( 14.530 . ( 14.630 we can verify ( 14.590 for g = 0. 

We finish the section remarking that the connection with this model will provide a 
useful tool in understanding the relation between Hurwitz numbers and Gromov-Witten 
invariants, also beyond genus zero, as presented in the Section 6. 



5 Critical behavior in the planar limit 

In this section we will show, by looking at the genus zero solution for the Gromov-Witten 
invariants obtained in the previous section, that there is a phase transition for p > 2 with 
the critical exponent typical of 2d gravity. 
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5.1 Critical properties of the planar free energy 

Since we have an exact expression for the genus zero Gromov-Witten invariants at all 
degrees, we can analyze the critical behavior by simply studying their asymptotic growth. 
We found, 

k\km{p-iy-i)k)V ^^-^^ 

up to a sign (— 1)^'^. By using Stirling's formula, we obtain 

No,k ~ e'^'^k-^/^, k^oo, (5.2) 

where tc is given in fl4.38p . Comparing with (12.31) we deduce that tc gives indeed the 
convergence radius of the expansion of the prepotential around t = oo, and we also 
deduce that 

7 = -^. (5.3) 

The above results are valid for p > 2. For p = 1,2 the series is convergent for alH > 0. 

We will now analyze the behavior of the prepotential near the critical point. By using 
the explicit relation between t and w (I4.39P we find 

w - Wc = A{e-*^ - e-^Y^^ + ■ ■ ■ (5.4) 

where Wc is given in (I4.37P and 

^-^^ (p_l)3-(P-ip p_i ■ (5.5) 

In order to extract the most singular part of the prepotential, we use (14.490 . which leads 
to the expansion 

= - log We - —{w -Wc) + — (e-*^ - e-*)i (5.6) 

ac^ Wc Wc 

Integrating this equation, we find 

which confirms the expected behavior (12. 4p and also gives us a precise value for the 
coefficient of the most singular piece. Again, the above holds only for p > 2, and for 
p = 1,2 the prepotential has a conifold-like behavior at t = 0. 

If we describe the saddle-point by the density p{h) given in (I4.43p . the critical behavior 
can be understood exactly as in the case of matrix models (see for example [35J). This was 
pointed out for the limiting case of Hurwitz theory in [21|, where F^^™^^^{fi) undergoes 
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the same kind of phase transition at /^c = e ^ For t > tc, the density p{h) behaves near 
the endpoint h = h as 

p{h) - I {b - h)^ , (5.8) 
where —7 = 1 — 6 and —f3 = 1 — a. Indeed, by expanding around this point we find 

00 

pih)-l = J2Mt,p){b-h)'--2, (5.9) 

k=l 

where 

If we now use the explicit expressions for the endpoints given in fl4.35p we find that, at 
the critical point, 

ai{tc,p)=0, a2{tc,p)j^0. (5.11) 

Therefore, criticality means, at the level of the density, that the leading branch cut sin- 
gularity near b is enhanced to 

p{h)-l^ {b-h)l. (5.12) 

The picture which emerges from this analysis is the following. In the large area phase 
t > tc, the planar model is described by the one-cut ansatz of the previous section, and 
the large radius expansion of the prepotential converges. At t = tc there is a phase 
transition controlled by the critical exponent 7 = —1/2, and not by the conventional one 
7 = which is found in the quintic, in local P^, and in other Calabi-Yau manifolds, as we 
reviewed in section 2. This is probably due to the fact that Zx^ is an equivariant partition 
function, and our result suggests that general equivariant partition functions like (13.121) 
will exhibit different critical behaviors for different choices of the equivariant parameters. 
On the other hand, in the small area phase t < tc the one-cut ansatz is no longer a valid 
solution of the planar model, since the equation f l4.34p does not have real solutions for w. 
This phase should be described by a chiral two-cut ansatz, as argued in [21] in the case of 
Hurwitz theory. It would be very interesting to obtain an explicit solution for this phase 
and to understand its geometric meaning. 

5.2 Multicritical behavior 

We have seen that, for p > 2, exhibits critical behavior with 7 = —1/2. In particular, 
if we regard this model as a deformation of the Hurwitz model, we have found that the 
critical behavior found in [21] at p 00 persists for all p > 2. 

It is then natural to ask if one can find multicritical behavior a la Kazakov [6IJ, and 
obtain the critical behavior of the (2, 2m — 1) models, with m > 3. The answer is yes, 
provided we turn on higher Casimir operators. This does not have a clear interpretation 
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in the context of topological string theory on toric Calabi-Yau manifolds, but it is natural 
to do if we view Zxj, as the partition function of 5d, M = 2 Abelian gauge theory (the 
5d version of [45j). We will however perform a precise analysis only in the undeformed 
model fl3.18p and at the planar level. It is clear that the deformed model will have the 
same multicritical behavior, but it is more difficult to solve in an explicit form. 

Let us then turn on higher Casimir operators in the Hurwitz model and consider the 
partition function 

^Hu.wit.(r.) = X: (^) 'iV^^(«)e--^(^)--^-/^^-^^^3 ^c.^ ^5_^3^ 

where 

N 

Cu = Y.hl (5.14) 

1=1 

Other definitions for the operators are possible. For example, using the operators 
defined in [IS] we obtain the partition function of topological string theory on with 
descendants of the Kahler class, and by using the operators defined in jlS] we obtain 
the partition function of f/(l), four-dimensional M = 2 supersymmetric gauge theory 
perturbed with Casimir operators. In any case, with the above modification, the saddle 
point equation becomes 

[ ^^'fri^ = + \og{h - 6), (5.15) 

W'{h) = T2{h-l) + n + Y,nh''~\ (5.16) 



where 



fe>3 



The resolvent is 



ujQ^h) = -W'{h) - -M{h) ^{h-a){h-b) + log(a - 6) + log /i - 2 log [Vh^ + y/h^ , 
^ ^ (5.17) 



where 



dz W'(l z) 1 
M{h) = (b — ^ ^ , , . 5.18 

To 2m 1-hz ^(1 -a^)(l -te) 



The density of eigenvalues is 



We can expand it around the endpoint h = b, and we find as before 

oo 

p{h) - 1 = Y,Mrk,p){b - hf--^. (5.20) 

k=l 
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The m-th critical point is achieved when 

p{h) - 1 {b - h)""-^ , m>2, 
i.e. we have to fine-tune the in such a way that 

Oik{rk,p) = 0, 1 < k <m-l. 
When this is the case, we get a critical model with 



(5.21) 



(5.22) 



^ m 



(5.23) 



To guarantee this, M{h) has to be a polynomial of degree at least m — 2 in h, which we 
denote by Pm{h), 

M{h)^Pm{h). (5.24) 



This is simply obtained by taking the first m — 2 powers of /i — & in the series 



Pmih) 



^ h-h\'^^(l 1 3 h-h 



We have used in the above that 



TT 



a — b 



fe=0 



2' 2' 2' a - 6 



(5.25) 



m-2 



— COS 



(2A:)! 



22fc(A;!)2(2A;+ !)■ 



(5.26) 



It is now easy to extract the derivative W^{h) of the m-th multicritical potential as 



w:,{h) = 4 p^h) 



h 



a — b 



1 - 



h 



1 - 



h 



(5.27) 



where the subscript + means as usual that we take the positive powers of h in the expan- 
sion of the above functions. We now list the first few critical polynomials. Denoting 



h-b 



(5.28) 



we have the simplified formula 



<(^) = 4 P^(^)^J 1 



(5.29) 
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For example, one has 




(5.30) 



We then see that, in what concerns critical behavior, and at least at the planar level, 
the general partition function ( ]5.13p is equivalent to the one matrix model, since one can 
reach all the (2,2m — 1) points. As we will see in the next section, the double-scaling 
limit of the critical model where the first and the second Casimirs have been turned on 
is the (2, 3) model (i.e. pure gravity). It is natural to conjecture that the double-scaling 
limit of fl5.13p is equivalent to the double-scaling limit of the one matrix model, and one 
can use it to describe the general (2, 2m — 1) models coupled to 2d gravity. This should 
also hold for the deformed model fl3.18p with higher Casimirs. 

On the other hand, Okounkov and Pandharipande have shown in [23] that the asymp- 
totics of the Hurwitz model (13.381) is equivalent to the asymptotics of the edge-of -the- 
spectrum matrix model, which in turn is equivalent to topological gravity in two dimen- 
sions. It is easy to see that the partition function (15.131) . when expanded in powers of 
the couplings r^, can be also expressed in terms of Hurwitz numbers. Indeed, it is closely 
related to (13.381) by a nontrivial map relating Tjt and Trl^^. In this sense, we can regard 
the Hurwitz model as an analog of the Gaussian matrix model with brane insertions, 
while the model (15.131) is an analog of the one matrix model. The relation between them 
and their equivalence in the double-scaling limit should be interpreted as some sort of 
open/closed duality for Hurwitz theory, in the spirit of [211 [25] . 

6 Higher genus analysis and double— scaling limit 
6.1 Explicit expressions at higher genus 

We have seen that the planar free energy of topological string theory on Xp undergoes a 
phase transition at small volume, with the same critical exponent as 2d gravity. We would 
like to understand now the behavior at higher genus. First we have to see if the Fg^it) 
have a critical point at the same tc-, and if the critical exponent depends on the genus as 
in (I2.12p . If this is the case, we can define a double-scaled theory at the transition point 
capturing the all-genus behavior, and we can try to extract the coefficients of the leading 
singularities at every genus in a unified way. 

We have then to compute Fg ^{t) from the partition function (13. 9p . Unfortunately, 
to the best of our knowledge, there is no systematic way to compute corrections to the 
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saddle-point from sums over partitions. We have then to use a different approach. We 
will in fact proceed in reverse: we will present an ansatz for the structure of Fg ^ (t) which 
manifestly has the critical behavior in (12.121) . Then, we will give evidence for the ansatz. 
The main evidence comes from the p —>■ oo limit. In this case, and as we explained in 
section 3, our model is equivalent to Hurwitz theory, and the ansatz holds thanks to the 
results of [22], [28], which were the source of inspiration for the conjectural expressions we 
will present. 

Our ansatz gives closed expressions for F-^ ^(t) in terms of a simple logarithmic ex- 
pression, 

F^^ = log(u; - ^c) - ^ log(p - 1) + ^(/ - 2p + 3) log w, (6.1) 
and for the higher genus Fg ^(t) as rational functions of the variable w 

^ '^■^ i=l 

The ttg^iij)) have the form 

'^M = ^1^, (6.3) 

where n is a positive integer and bg i{p) are polynomials in p with rational coefficients. We 
have verified this ansatz by direct computation up to genus 4. In this ansatz the Fg ''(t) 
are determined by 5((7 — 1) coefficients ag^i{p), which can be uniquely determined from the 
genus g Gromov-Witten invariants up to degree 5{g — 1). These can be computed from 
the explicit expression (13. 9p for Zxp- One then verifies that the expression (16. 2p obtained 
in this way reproduces correctly the Gromov-Witten invariants of higher degree. 
For genus 2, for example, one finds in this way: 

1 Pip -2) 



02,5 (P) 



2880 (p- 1)2' 

1 12-14j9 + 7p2 
~2880 (p - 1)4 ' 



1 36 - 106p+161p2- 204p3 + i7ip4_72p5 + 12/ 



«2,2(j3) 
0'2,l{p) 



1 36 - 90p+121p^ -60p^ -bp" + 12p^ -2p' 

2880 (p-iy^ 

1 1 



240 (p- 1)10 

Once the coefficients ag^ip) have been obtained, one can deduce closed formulae for the 
Gromov-Witten invariants Ng^d for all d by using Lagrange inversion. We give the results 
for g = 1 m Appendix A. 
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It is interesting to notice that the ansatz fl6.2p is very similar to the holomorphic 
ambiguity in the standard B-model topological string [3]. There, one finds that Fg{t) is 
given by a piece which is completely determined in a recursive way by special geometry 
data and amplitudes at lower genera, plus an undetermined, holomorphic piece of the 
form 

M^) = (6-5) 

where z is a natural coordinate in the complex structure moduli space, A{z) is the dis- 
criminant locus, and Pg{z) is a polynomial in z with unknown coefficients. The exponent 
2 — 2(7 in the discriminant is related to the critical exponent 7 = characteristic of the 
c = 1 behavior fl2.18p . In our model, the F^''(t) is given entirely by an analog of the 
holomorphic ambiguity, with a different singularity due to the different critical exponent. 
This also confirms the role of w as a natural mirror coordinate, and that the relation 
between w and e~* is a mirror map. 

To further justify this ansatz, let us first recall the results of [221, [28] for the generating 
functionals of simple Hurwtiz numbers: 



F,"^(/x) = -l(log(l-x)+x), 



P (y) ^^-^^ 



1=1 



Here, x is the variable introduced in (I4.60p . which is related to fi by ( 14.62^ . Moreover, 
there are explicit expressions for Cg^i in terms of Hodge integrals. In particular, one has 
that 

which is a correlation function at genus g in 2d topological gravity (see for example [35] 
for a review). This result will be useful later in this section. 

It is easy to see that, in the hmit fl37[5|) - fl3A6| ). (Kl^ becomes and that 

the ansatz (16.21) leads to expressions which are compatible with (16.61) . provided the poly- 
nomials V{w,p) satisfy certain conditions. These conditions can be written as limiting 
conditions on the coefficients ag^i{p): 

hm /(^'-i)"2's.,(p) = i-iycg,i. (6.8) 

p— >oo 



Using (14.651) as well as the result 



w -Wc^ (6.9) 
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in the limit fl3.15l) - fl3.16l) . it is easy to see that, if (16.81) holds, then fl4.59p is obtained. For 
example, by using the above expressions for genus two, one obtains in the p oo limit, 

1 X . ^ 1 / 12 o 2 A , , 

^F^^it) + X , 6.10 

p2 2 V ; (1 -;^)5 \^2880-^ 2880-^7' ^ ' 

which is the expression found in [221 [28] . 

6.2 Double— scaling limit and Painleve I 

Using the above results for higher genus, it is immediate to see that the ansatz (16.21) 
gives the critical behavior we are looking for. This is a consequence of f l5.4p . We want to 
analyze now the coefficients of the leading singularity. The behavior of the singular part 
of the planar free energy (15.71) suggests defining a scaled string coupling z as 

z'/^ = 97'^^^{w - = 9;'^^^A\e-^^ - e'^i (6.11) 

We can now consider the double-scaled theory in which we take the limit 

t tc, — > 0, z fixed. (6-12) 
In this limit, the total free energy of the model becomes the double-scaled free energy 

Fx,^F^{z). (6.13) 

Taking into account (15.71) and the expression (EH) for "(t), we find that up to genus 
one 

Fds(z) = -—z^'^ - — log z + ■ ■ ■ . (6.14) 

dsv ^ 15 48 ^ ^ ^ 

This is, up to this order, the perturbative expansion of the free energy of 2d gravity, 
-^(2,3) (-2^)- We recall (see [33]) that -F(2,3)(-2) is determined as a function of z by the following 
equation, 

^{2,3)(^) = -^(^), (6.15) 
where u{z), the specific heat, is a solution of the Painleve I equation 

- ^m" = z (6.16) 

with the asymptotics at weak string coupling 

u{z) = 2;2 + ■ ■ ■ , z ^ oo. (6-17) 
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^^^^-^) = ^7;:3w ^^2. (6.21) 



This leads to the asymptotic expansion at large 2, 

Fi2M^) = -U'^ - ^ log - + E (6.18) 

g>2 

with 

7 245 259553 

ao = , 03 = , 04 = , (0.19) 

5560 331776 159252480 ^ ' 

an so on. 

In view of the above results for genus = 0, 1, it is natural to conjecture that the 
double-scaled free energy of topological string theory on Xp equals the free energy of 2d 
gravity, at least in the genus expansion 

i^ds(^) =i^{2,3)(^). (6.20) 

This means that the coefficient of the most singular part of Fg ^(t) is identical to the 
coefficient ag in f l6.18p . when expressed in terms of the variable z defined in (16.111) . In 
terms of the ansatz (16.21) . in order to verify the conjecture (I6.20p we have to verify that 
the polynomials Vg{w,p) have the following value at the critical point w = Wc 

It is easy to check with (16. 4p that this is indeed the case for g = 2, and we have verified 
(16.211) up to genus four. We can indeed provide evidence for (I6.2ip at all genus in the 
limit p ^ 00. Taking into account the overall factors of p, one easily finds that in this 
limit the conjectured equality (I6.2ip becomes 

Pg{l) = A^'^ag, g>2, (6.22) 

where Pg{x) are the polynomials that appear in the Goulden-Jackson- Vakil expressions 
for ^7^^™*^ in (16. 6p . But due to (16. 7p this reduces to checking the following equality for 
2d topological gravity correlators, 

^^^^(ar'). = 4^~'%, 9>2, (6.23) 

which indeed was shown to be true in ^62j. This verifies our conjecture in the limiting 
case j9 — i> cx) at all genus. 

To conclude this section, we have provided strong evidence that the critical behavior 
of topological string theory on Xp at small radius is controlled by the Painleve 1 equation 
for all p > 2 (the dependence on p only enters in the normalization factor relating z to 
the "bare" couplings t,gs). Therefore the double-scaled theory at the transition is 2d 
gravity (i.e. the (2,3) model). This is in contrast to conventional topological strings, 
which undergo a phase transition at the conifold point controlled by c = 1 string theory 
at the self-dual radius. The possibility that different kinds of singularities in Calabi-Yau 
manifolds lead to different universality classes was pointed out at the end of [1]. Here we 
have found a concrete realization of this possibility. 
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7 Non— perturbative proposal and 2d gravity 



7.1 The non— perturbative proposal and its critical properties 

Topological string theories, as well as ordinary string theories, are defined in principle 
only perturbatively, in a genus by genus expansion, and it is natural to ask whether 
a nonperturbative definition is possible. In the case of the backgrounds that we are 
considering, i.e. Calabi-Yau manifolds of the type Xp, it was proposed in ^ based on 
previous works [63l[67] that such a definition can be given in terms of a q-deformed version 
of two-dimensional Yang-Mills theory with gauge group U{N). The strong coupling 
expansion of this partition function is given by the expression 

Z,YM = $^(dim,i?)^g^^^(«)/2e^^^(^). (7.1) 

R 

Here, R is a. tableau representing an irreducible representation of U{N), C2{R) = k,r + 
Ni{R), and the quantum dimension of R is given by 

l<i<j<N LJ J 

In [H] it was shown that the asymptotic expansion of (17. ip can be written as 

ZqYM = ^Y1 Zr,R2^^ + P9sl)ZR^R^{t - pQsl). (7.3) 

In this expression, t is given in terms of the gauge theory data as 

t=]^{p-2)Ngs-ie. (7.4) 

The amplitude Z^^j^^ is a perturbative open string amplitude on the Xp Calabi-Yau and 
in the presence of two stacks of D-branes, labeled by the representations i?i,i?2- An 
explicit expression for Z^^j^^, written in terms of the topological vertex, can be found in 
[6]. When Ri = R2 = • is the trivial representation, one recovers the perturbative closed 
string amplitude 

Zt, = Zx,, (7.5) 

up to minor overall factors that do not affect the worldsheet instanton expansion. There- 
fore, the nonperturbative completion (17.11) gives the square of the perturbative answer (as 
expected from [63]), plus some extra insertions of branes. 

Later on, it was noticed in [T2| IT3| [H] that, at large N, the partition function (17.11) 
(or more precisely, its zero charge sector I = 0) has a third-order phase transition for 
p > 2 which is qualitatively similar to the Douglas-Kazakov transition of 2d Yang-Mills 
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theory. In terms of the parameter t which is identified with the Kahler parameter of the 
perturbative theory, and for ^ = 0, the critical point is given by 

tnAp) = \p{P - 2) log(^l + tan2(^)) • (7.6) 

Based on the analysis of [13], we expect that this critical value gives the radius of con- 
vergence of the strong coupling expansion of (17. ip . in the same way that (I4.38P gives the 
radius of convergence of the perturbative theory. 




5 10 15 20 25 30 



Figure 3: The curve in the top shows the critical value of the Kahler parameter in the 
perturbative theory, which goes as 2 log p for large p. The curve in the bottom shows the 
critical value in the q-deformed 2d Yang-Mills theory, which asymptotes 7r^/2 as p goes 
to infinity. 

Therefore, hothihe perturbative theory, defined by Zx^,-, and the nonpertubative defini- 
tion proposed in [6], undergo a phase transition at small radius for p > 2. It is interesting 
to compare the critical behaviors as a further probe of the proposal of E] . The first 
things to compare are the radii of convergence tc{p) and tnp(p), as a function of p. In 
Fig. [3] we plot the curves (14.381) and (17. 6p for 2 < p < 30. In the perturbative theory, one 
has 

tc 2 log p, p-^ oo, (7.7) 

while 

^np^y, p^oo. (7.8) 

Therefore, at large p, the proposed nonperturbative completion has better convergence 
properties function of t at small radius. 

Another interesting difference between both models is the universality class of the 
theory at the transition. We have seen in the previous section that one can define a 
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double-scaled theory at the critical point, and further analysis shows that this theory is 
2d gravity and it is governed by the Painleve I equation. A similar study in the case 
of q-deformed 2d Yang-Mills theory is still lacking, but the undeformed theory (which 
appears as p ^ oo, as explained in [121 [131 E]) has been analyzed from that point of 
view in [18j. It turns out that there is a well-defined double-scaled theory as well, but 
this time it is in the universality class of the Gross- Witten-Wadia unitary matrix model 
[151 [IB]- The double-scaled free energy F{z) is given by F"{z) = v{zY/4:, where v{z) 
satisfies the Painleve II equation [T7] . 

2v" -v^ + zv = 0. (7.9) 

As argued in [19j, this also describes the universality class of pure 2d supergravity. As in 
the case of the perturbative theory, where the universality class is the same for all p > 2, 
we expect that q-deformed 2d YM has a double-scaled theory described by Painleve II 
for all p > 2. Some evidence for this was given in [T3], based on the behavior of the 
instanton suppression factor. It was shown by Crnkovic, Douglas and Moore ^20j that 
the Painleve II equation (17. 9p has a unique real, pole-free solution which gives a non- 
perturbative solution of the double-scaled Gross-Witten-Wadia unitary matrix model 
and of 2d supergravity [19J. This solution, which we will call vcbm, is also the correct 
nonperturbative description of the double-scaled theory at the critical point of q-deformed 
2d Yang-Mills. 

We therefore reach the conclusion that perturbative topological string theory on Xp 
exhibits a phase transition in the universality class of 2d gravity, while the proposed non- 
perturbative completion exhibits a phase transition in the universality class of 2d super- 
gravity. In the remaining of this section, we will present some comments and speculations 
on the possible implications of this fact for 2d gravity theorie^. 

7.2 Implications for 2d gravity 

It is well-known (see again |i35j) that the Painleve I equation does not define 2d gravity 
beyond the perturbation regime, since the resulting series for the specific heat is not Borel 
summable. There have been various proposals to provide a nonperturbative definition of 
2d gravity which agrees with the perturbative expansion (I6.18P at large z. A particularly 
interesting proposal was made by Dalley, Johnson and Morris in [M]. By requiring that 
the KdV flows of 2d gravity hold nonperturbatively they found a modified string equation, 

uTZ^ -^rnz" + ^{n'y = 0, n = u^-u"/3-z. (7.10) 

Notice that 7^ = is just the Painleve I equation, with a different normalization than 
(16.161) . In this normalization, the free energy computed from u is actually twice the free 

^ We would like to thank Clifford Johnson for pointing out to us the relevance of [SU [65] to our work, 
and for an illuminating email exchange. He informed us that he had previously envisaged the possibility 
that the results of [Ml [65] might be relevant in topological string theory. 
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energy of pure gravity (see for example [35] )• It can be shown that the above equation has 
a smooth, real solution u^jm, which at z ^ cxd reproduces the genus expansion obtained 
from Painleve I, and vanishes as l/z^ when z — ^ — oo. This solution was proposed in [61] 
as a nonperturbative definition of 2d gravity. 

The equation (17.101) can be generalized by including an extra parameter F, as follows: 

uTZ^ - -nn" + -in'? = r^. (7.11) 

2 4 

It was argued in [65] that this parameter corresponds to the introduction of an open string 
sector. This equation has solutions that are open string generalizations of mdjm and which 
we denote by Mdjm- For ^ ~^ ~oo they behave as 



r2 _ 1 

«LmW ^- (7.12) 



It was shown in [65] that any solution v to the Painleve II equation (17. 9p can be mapped 
to a solution u of (17. lip with F = ±1/2, through the so-called Miura transformation (see 
also [66]). In particular, the solution vcdm maps under the Miura transformation to 
the solutions m^jm^''^- Moreover, if we denote by Zpn and the partition functions 
associated to the solutions related by Miura transformation, we have the relation 



jpii = ^r=i-^r=-i- (7-13) 



Notice that this is very similar to (17. 3p . and in both cases one needs to introduce open 
string sectors (represented in (17.31) by the D-branes and here by the F parameter). This 
seems to indicate that the relation ( I7.3P between the perturbative topological string am- 
plitude and its nonperturbative completion is inherited as the relation (17.131) between 
their critical counterparts. 

Recall that the proposal of [HI |63l |67] is to regard q-deformed Yang-Mills as a nonper- 
turbative completion of topological string theory on X^. If we focus on the critical points 
of the corresponding theories, the proposal indicates that the nonperturbative completion 
of the 2d gravity perturbative expansion (16.141) should be associated to the solution f cdm 
of Painleve II. In other words, the connection to Painleve II should provide the extra non- 
perturbative information needed in 2d gravity. As we have just reviewed, wcdm is related 
to the solutions u^^^"^ with an open string sector, as encoded in (I7.13p . However, the 
perturbative theory is obtained properly speaking by "removing" the open sector, as in 
(17. 5p . We conclude that the embedding of 2d gravity in perturbative topological string 
theory on Xp, together with the proposal of [HI [HSl EZ], suggest that the natural nonper- 
turbative completion of 2d gravity is indeed the solution presented in [61] itojM = ^djm- 

We emphasize that, although we find this line of argumentation suggestive and com- 
pelling, it should be put on a much firmer ground by a better understanding of the 
nonperturbative proposal of [6l [63l |67] and of the embedding of the critical theories. 
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8 An open— closed duality 



In section 3 we saw that the partition function Zxj, is related to a configuration involving 
a topological D-brane in with framing f = p — 1. In this section, we will point 
out a different relation which holds at genus zero and generalizes an open-closed duality 
proposed in {,29| in the context of M-theory. 

Consider then the open string background given by a D6 brane wrapping a Lagrangian 
submanifold with topology x in and with trivial framing / = 0. As explained 
in |6S1 ESj, this leads to a superpotential W in four dimensions which is given by 

d=i 

where u is an open string modulus. From the point of view of the open string amplitudes 
Fgj: which appeared in section 3, this superpotential is the generating functional of disk 

amplitudes with g = and a vector k with one single entry k^ = 1. As pointed out in 
[29] . this superpotential is related to the free energy of the resolved conifold by 

Wf=o (8.2) 

after identifying t = —u. This relation between superpotential and prepotential is typical 
from open/closed dualities, where the presence of the brane in the open side is traded by 
the presence of flux in the closed side. 

Consider now the same brane configuration but with arbitrary framing /. The super- 
potential has been computed in ^69j and reads 

Wf = Y.Wf^,e'-, Wf,u=^-^^\{{kf-3). (8.3) 
fc=i ■ j=i 

Comparing this to our prepotential Fq ^{t) given in fl4.53p . we find 

^ -'^^ (8.4, 

with the choice of framing 

f = {p-lf. (8.5) 
In particular, the genus zero Gopakumar-Vafa invariants riQ^dip) are given by 

noAv) = \NdU), (8.6) 

where Nd{f) are the open BPS invariants associated to disk instantons. Notice that 
this relation is very different from ( 13.46^ . The invariants Nd{f) can be written as linear 
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combinations of BPS invariants N{i^g=o{f) involving hook tableaux of d boxes. But the 
open string backgrounds involved in fl3.46p and in (18. 4p are different, since in the first case 
one has framing f = p — 1 while in the second case one has / = {p — 1)^. 

The result f l8.4p indicates the existence of an open-closed duality between closed string 
theory on Xp and open string theory on in the presence of a framed D6 brane with 
framing (18.50 . generalizing in this way the relation (18.20 for the resolved conifold p = 1. 
In the p = 1 case, the open/closed duahty can be derived by lifting both configurations to 
M-theory on a manifold with G2 holonomy and with topology §^ x [29]. For general 
p, one should be able to lift the D6 brane with framing to M-theory and relate it there to 
closed string theory on Xp, explaining in this way the equality (18.41) . This would produce 
a new and interesting class of M-theory dualities. 

9 Conclusions and open problems 

The main conclusions of our analysis are the following: 

• We have shown that, in some backgrounds, topological string theory exhibits a 
critical behavior which is not in the universality class of the c = 1 string, and leads 
to a double-scaled theory which is identical to pure 2d gravity (in perturbation 
theory) . 

• From a more methodological point of view, we have shown that double-scaling limits 
can be useful to characterize critical string theories, and they allow us to partially 
capture their all-genus behavior. This is a useful strategy to address questions 
related to their nonperturbative behavior, and we have argued in this paper that 
the double-scaled theory reflects aspects of the nonperturbative structure of the 
original theory. In particular, we have seen that the proposal for a nonperturbative 
completion of topological string theory on Xp leads to a double-scaled theory with a 
well-defined nonperturbative description in terms of the solution fcDM to Painleve 
II. 

• Conversely, we suggested that this result gives further support to the proposal of 
[M] for a nonperturbative completion of 2d gravity. 

Our work also leaves many open questions, and we list some of them here: 

• From a technical point of view, our analysis has some gaps which should be filled. 
For example, it would be interesting to give a first principles derivation of the 
higher genus ansatz (16. 2p . as well as a constructive way of computing the coefficients 
ag^i{p). In order to do that, one should develop general techniques to compute higher 
genus corrections for sums over partitions, which would find applications in many 
other problems related to topological strings and instanton counting. A systematic 
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method to treat higher genus corrections would also make possible to derive the 
precise properties of the double-scaled theory, as in the case of matrix models. 



We have not determined the nature of the theory in the phase with t < tc- In 
principle, this corresponds to a two-cut solution of the saddle-point equations, but 
finding it is technically challenging. In this respect, it would be very interesting 
to develop mirror symmetry techniques to analyze this model (both at the planar 
level and at higher genus). Some preliminary steps in this direction have already 
been given in [IQ], but much remains to be done. This would also shed light on the 
B-model analysis of equivariant topological strings introduced in 



• We have found evidence for an open-closed duality involving a lift to M-theory, 
and generalizing the observations of [29]. It would be very interesting to use G2 
holonomy manifolds and justify the relation between the D6 brane background with 
framing / = (p — 1)^ and the closed Xp background. 

• There have been proposals for nonperturbative formulations of topological string 
theory on other toric manifolds, like local |70]. It would be very interesting to 
analyze the critical behavior of these models and their universality class. Based on 
the results of this paper, we would expect a "well-behaved" critical theory, probably 
related to doubly-scaled unitary matrix models. 

Note added: After this paper was submitted, the work [71j appeared, which studies 
the B-model realization of topological string theory on Xp. They derive the following 
explicit expression for the mirror map. 

It can be easily seen that after setting (—1)^ = w — 1 the above expression becomes 
precisely our equation fl4.34p . This confirms our conjecture that the relation between w 
and t is indeed the mirror map of this model. By using B-model techniques, [71] also 
derives an expression for Fi which agrees precisely with (16.11) . 
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A Lagrange inversion 

The Lagrange inversion method makes possible to exphcitly invert a relation between two 
variables in power series form. A good exposition can be found for example in \72\ . 

Let z = f{w) be a relation that implicitly defines w as a function of z around w = 0. 
We can assume that /(O) = 0. We can invert this relation and find the inverse function 

w = r\z) (A.2) 

as follows 



oo 

W 



V-Resf— VV", (A.3) 



where the residue is computed around w = 0. A slight modification of this result makes 
possible to compute the function 

g{w) = g{f-\z)) (A.4) 



as 



n=l ^ 



We assumed again that g{0) =0. 

Let us now apply this procedure to the equation (14.340 relating t and the variable w 

g = w(P-l)'(u;-l-l), (A.6) 

where q = e~*. Set f = {p — 1^ , x = w — 1. The equation defining x is 

q = -x{l + xy-\ (A.7) 
We want to find now x (or w = 1 + x) as a function of q. By flA.3l) . we have to compute 

n (l + a;)"(/-i) n ^ k\{n{f - 1) - 1)\ ^ ' ' v ; 

which equals 



k=2 



Therefore, 

oo n—1 



W 

n=l k=l 

where the n = 1 term is taken to be 1. 



oo ft— 1 / ^ x„ 
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Lagrange inversion can be also used to find explicit expressions for the Gromov-Witten 
invariants of Xp from (16.11) and fl6.2p . Let us calculate the genus one invariants. In terms 
of X = w — 1, we have 



Fi = loe 

24 ^ 



1 + (p - ifx 



+ ^(/-2p + 3)log(l + x). 



Setting again f = {p ~ 1)^ we have to compute the residue: 



-1)" 



Res 



2An 

and we easily obtain 

^=-1 ck-£ 



f 



f + 2 



X" 1 + /X (1 + X"' (1 + 



'v.. = 4 E ^ n('=(/ - 1) + J- - 1) - Y Jm-m <^f + 



24A;^ 

£=0 j=l 



We can check the limiting behavior fl3.34p 



lim 

p^oo p^^: 24/c 



^ k-2 

-T 

Ah 



i=0 



(A.ll) 



(A.12) 



(A.13) 



(A. 14) 



which is indeed the expected answer. 



B Useful Integrals 

In this appendix we shall summarize some integrals that have been used to perform the 
large N analysis in sec. 4. 

All the the integrals appearing in the computation of the resolvent uj{z) can be reduced 
to the following elementary indefinite integral 

' (B.15) 



w 



1 J / (^^y{w-e-f^)(s-e-f}+^/{s-e-l^){w-e-'') 

The computation of the planar free energy involves instead the following families of definite 
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integrals. First we have to consider 
d(f) log 



2 



2,e^ + e"' — (^e'^ — e^) cos ( 



1 r 

2 7-. 



d(f) 



(e-/3/2 + e-^/2)' 
log ^ ^ + log 1 - 



e/3/2 _ e7/2 
e/3/2 + e7/2 



e I + 



(B.16) 



+ log 1 - 



e/3/2 _ e7/2 
e/3/2 + e^/2' 



= 47rlog^ 



-/3/2 



+ e- 



-7/2 



+ 27rLi2 



e/3/2 _ e7/2 
e/3/2 + e7/2 



This integral is actually computed by expanding the square, then by expanding in series 
the logarithms and finally integrating the series. The final series can be summed in terms 
of polylogarithms and logarithms . 
With the same technique, we can also show that 



^2^\f d0 V (1 - Ce'^) log (1 - Ce-^'t') = 



= -TT (log (C') log' (1 - C') + 2Li2 (1 - C') log (1 - C') - 2Li3 (l - C') + 2C(3)) 

(B.17) 



and 



2 

+ 67rLi2 



log^ {A -Bcoscj)) = TTlog^ ( ^ + ^VA^-B^ ] + 



D 



A + VA^ - 



) jlog (^^ + lVA2-B2^ +6J2 



D 



A + VA^-B\ 

(B.18) 



Using the result for I^, we can obtain 
1 



1 d r 

X4 = - / d(p log' {A-B COS0) /{A - B COS0) = -— / d0 log^ {A-B cos0) 
2 6 oA J _^ 



TT log 



2 f 2(b2_a(aWA^-B^))' 



+ 2Li2 



(B.19) 

which can be used to evaluate the following integral appearing in our computation of the 
planar free energy 

r , ^2 ( + e"^ - (e-^ - e-^) cos ( 

Jo 



d(j) 



((e-^ - e-T) COS0 + e-/^ + e^^) 



log' 



= TT 2 log' 



'(e/3/2 + e^/')' 



' ' ' " ^ ' + Li2 ( tanh' ' ^ 



(B.20) 
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The results of these integrals can be also presented in an equivalent and more useful way, 
if we use the identities 



7r2 



log(l - x) \og{x) + Li2(l - x) + Li2(x) - — = 0. (B.21) 





and 



11 TT^ fw — 1\ 

LisH = -\og\w)--\og{l-w)\og\w) + —\og{w)-U^{l-w)-U^ +C(3). 

5 2 5 \ w / 

(B.22) 
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